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ABSTRACT 

We propose a simple explanation for the constant mean brightness observed between out- 
bursts in dwarf novae. Secular changes in the total energy dissipation rate of the accretion 
disc brought about by variations in surface density, temperature and disc radius can be reg- 
ulated by the gradual cooling of a small, critically-stable hot inner region. The hypothesis is 
supported with two-dimensional time-dependent numerical models of dwarf nova accretion 
discs. 

Key words: accretion, accretion discs - instabilities - binaries: close - methods: numerical - 
novae, cataclysmic variables. 



1 INTRODUCTION 

Dwarf novae have been observed extensively since the discovery of 
U Geminorum in the mid-nineteenth century. They a ccount for well 
over half of the non-magnetic cataclysmic variables iRitter & Kolbl 
t_1998), and consist of a white dwarf accreting from a low-mass star 
via Roche lobe overflow. The dwarf novae are observed in two 
states: a quiescent, almost constant-luminosity low state that may 
last for months, and a regular (but aperiodic) high luminosity, out- 
burst state lasting for several days. The transition between the states 
is fairly rapid; typically the V magnitude changes by four magni- 
tudes over the course of a day or so. 

The outbursts are attributed to a thermal-viscous instabil- 
ity in the accretion disc, often referred to as the disc instability, 
which is associated with the steep dependence of opacity with 
temperature at the point where hydrogen begins to ionize (Th ~ 
6500 K). Our understanding of the physical properties of the ac- 
cretion disc in quiescence is less complete. In particular, while 
great strides have been made in the identification of a magneto- 
rotational in stability (MRI) responsib le for the transport of angular 
momentum I'Balbus & Hawlev"l99l'), the efficiency of this insta- 
bility is se nsitive to the degree of ionization in the accretion disc 
jMatsumoto & TaiimaM199 5). It remains unclear whether the MRI 
can sustain the necessary angular momentum transport appropriate 
to a cool, quiescent disc. It is po ssible that oth er viscosity mecha- 
nisms are at work in quiescence jMenoj2000l) . 

There have been several rigorous one-dimensional calcu- 
lations of dwarf nova outbursts wi t hin the framework of th e 
disc instability model ^Smak" Il984bl) : llchikawa & Osakil h992h : 
[Hameuryetal. ( 1998) and many others). The recent review of 
iLasotal J200ll) describes the current state of the disc instability 
model, and contains an excellent commentary on its scope and limi- 
tations. In his review, Lasota describes quiescence as 'the Achilles' 
heel' of the model, because the model handles the behaviour of 



dwarf novae in quiescence rather poorly. One of the major prob- 
lems with the numerical models has been that when the correct 
boundary conditions are used (in particular an outer disc radius that 
is not constrained to be fixed) the models are unable to reproduce 
the observed constant level of brightness during quiescence. The 
brightness almost always increases between the outbursts. 

Visual observations of dwarf novae by amateur astronomers 
cover more than a century in some cases and decades in many oth- 
ers (for example ICannizzo & MatteH il992h ). These visual obser- 
vations show scatter at the level of 0.5 to 1 magnitudes, so it is 
difficult to see short-term variations, but it is clear that the mean 
visible flux remains approximately constant throughout quiescence. 
In contrast, disc instability models consistently predict an increas- 
ing brightness during quiescence. The amplitude of the predicted 
increase lies in the range 1 to 3 magnitudes, and would certainly 
have been observed if present. 

In this paper we suggest an explanation for this discrepancy. 
In the next section we present an analytic model of how the vis- 
cous energy generation rate responds to variations in the physical 
state of the disc over the quiescent interval in a dwarf nova. We 
show that in a disc with a constant Shakura-Sunyaev viscosity, the 
total viscous energy dissipation rate increases in response to the 
expected changes in surface density, temperature and outer disc ra- 
dius. However, we show that it is possible for the total viscous en- 
ergy generation rate to remain constant if there is a small portion of 
the disc that remains in the high-viscosity state during quiescence. 
This hypothesis is supported by two-dimensional smoothed parti- 
cle hydrodynamics (SPH) simulations of accretion discs in Section 
3, which incorporate the evolution of the free outer disc boundary 
in response to the full two-dimensional tidal potential of the binary. 
We discuss the possible implications for the disc instability model 
in Section 4. 
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2 QUIESCENT LUMINOSITY 

We consider the local viscous energy dissipation rate in an accre- 
tion disc, which is given by 



D(r,t) = -Ki)S(r,t)-^^^^3, 



(1) 



This definition differs by a factor of two from that which appears 
in some treatments. Here we express the local dissipation rate inte- 
grated over the vertical interval — oo ^ 2 ^ oo, not just from the 
mid-plane. We assume that the viscosity c an be described by the a 
prescription of lShakura 8i SunvaevH 19731) : 



u — ac^H = ac. 



GMi 



(2) 



where the viscosity parameter a is a dimensionless constant. 

Now, at a time t, the total luminosity due to viscous energy 
dissipation in the accretion disc is 

L{t) = / D(r,t) 27rrdr, (3) 
where i?in and Rd are the inner and outer radii of the disc. Then, 



L(t) cx / aT{r,t)T.{r,t)r{t)~^dr, 

J Rir, 



(4) 



where T{r, t) is the local mid-plane temperature at time t. 

Since the conditions in the accretion disc are not steady dur- 
ing quiescence, there is no exact analytic expression for the time- 
dependent evolution of the surface density, E(r-, t) and the mid- 
plane temperature, T{r, t). However, we do know that these quan- 
tities must increase during quiescence if another outburst is to fol- 
low via disc instability. Indeed, given that the critical surface den- 
sities for triggering a nd suppressing t he disc instability vary almost 
linearly with radius iCannizzo et alJi.1988) . one could envisage a 
model in which the surface density profile remains linear with ra- 
dius while its gradient or mean level steadily increases. Further, we 
could assume that the radial temperature profile remains almost flat 
during quie scence, as suggest ed by eclipse mapping of the dwarf 
nova Z Cha iWoodetalJigS^) . We do not make these assumptions, 
but for simplicity we express the product 



T(r,t)E(r,t)=AWr" 



(5) 



where A is independent of radius and n can take any value. Under 
the two assumptions discussed above, n would be unity, but we 
keep it as a free parameter to take into account other possibilities. 
Its exact value (positive or negative) does not affect our argument 
here, but we do assume that since temperature and surface density 
increase during quiescence, the quantity AA/At > 0. In this case. 



L{t) oc oLcA(t) 



R" 



(6) 



where the entire disc is on the cool branch of the disc instability 
curve, with a viscosity parameter etc- If we consider how the lumi- 
nosity will vary during quiescence, assuming that the inner radius 
remains constant. 
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The disc radius is expected to decrease slowly between dwarf 
nova outbursts, since the accretion rate onto the primary is small 



and material with low specific angular momentum is accumu- 
lating near the outer edge from the gas stream. This ha s been 
confirmed by one-dimensional disc instability calculations iSmakI 
J984a; Ichikawa & Osaki 1992). However, this term is dominated 
by the increasing surface density and temperature, and the total disc 
luminosity increases between outbursts. 

Now consider a disc with two regions: a hot, inner region on 
the upper branch of the instability curve with ah, and a cool outer 
region as before. We term the transition radius between the regions 
Rt. In this case, 

fRt 1 rRi , 

L{t)(x I QhAh(t)r"-2dr+ / acA,{t)r"'^dr, (8) 
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The second and fourth terms of this equation are identical to the 
two terms in equation Q however, there are now two additional 
terms involving the rate of change in Rt and Ay^. The presence of 
a hot inner region means that there are more parameters available 
to change dL/dt. As an example, the coefficients of the two new 
terms in equation [Tol involve the hot-state viscosity, Qh, while the 
coefficients of the other terms involve only Qc. Since we expect 
ah ~ lOctc in a dwarf nova accretion disc, any changes in dissipa- 
tion rate brought about by the evolution of the cool part of the disc 
could be offset by much smaller changes in the hot, inner region. 

We stress that the above analysis is in no way meant to be a 
self-consistent solution of the accretion flow. We simply demon- 
strate that a hot, high-viscosity inner region provides many more 
avenues for evolution of the disc. For example, the situation is com- 
plicated by the fact that we do not know the correct value of n or 
the correct form of A{t). In particular, it is extremely difficult to 
predict the response of Rd and Rt in a non-steady accretion disc. It 
is only in a numerical calculation that these effects can be studied 
in more detail. 

In the next section, we demonstrate the effect of maintaining 
a hot, high-viscosity inner region with two-dimensional smoothed 
particle hydrodynamics calculations of accretion discs in dwarf no- 



3 TWO-DIMENSIONAL SPH CALCULATIONS 
3.1 Numerical method 

We use a two-dimensional (2D) smoothed particle hydrodynamics 
scheme to model the evolution of an accretion disc in a close binary. 
The model includes the full tidal potential due to the secondary star 
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and has been discussed in d etail and tested with reference to dwarf 
novae m iTruss etalji2000l) . 

SPH uses an ensemble of particles to describe a fluid. The ad- 
vantages of such a Lagrangian scheme to model these systems im- 
mediately becomes apparent: the response of the free outer bound- 
ary of the accretion disc to the mass stream and the tidal potential 
are treated in a self-consistent manner. F or a wider d iscussion of 
SPH, the reader is referred to the review of Monaahan (1992). 

We solve the SPH momentum equation in the Lagrangian (co- 
moving) frame, with an artificial viscosity parameter ( for each 
particle, i, such that 



IT 



+ 



(11) 



where pij and Cij are the average density and sound speed of a pair 
of particles i and j respectively, Wij is an interpolating kernel (here 
taken to be a cubic spline) and 



Hvii 



(12) 



The softening parameter r) avoids singularities and H is the length 
scale over which the viscous energy is dissipated. Here, this is set 
to the scale height, H = We do not use a viscosity term that 
is quadratic in /iij , often denoted by the parameter /3 in SPH equa- 
tions. A linear term is quite sufficient in this case. .Murray 1.199^ 
has shown that in 2D, the C, term introduces a shear viscosity 



(13) 



so the implementation is equivalent to a Shakura-Sunyaev viscosity 
of a = ("/8. We also solve the thermal energy equation 



1 

2 ^ 



2 ' 1 



(pi. 



Pij 



(v.-v,).V,mj. (14) 



The thermal-viscous instability is incorporated using a viscos- 
ity switch for each particle based on the local surface density. The 
critical conditions for the tri ggering or suppression of the instability 
are roughly linear in radius iCannizzo et alll988l) . Here, we take 

'^max "^hot i 



r ^ 

^cold 1 

a 



(15) 



where Chot and Ccoid are constants and a is the binary separation. 
If the surface density of a particle exceeds Emax, then its viscosity 
is increased smoothly from a value <^c to a value i^h on a 'ther- 
mal' time-scale rtiiggcr, where rtiiggcr << tviscous. Its viscosity 
remains high until the local surface density falls below Smin, when 
it is returned to C,c on a similar time-scale. In contrast to previous 
implementations of this code which were fully isothermal, here we 
change the sound speed of the gas whenever the critical surface 
density thresholds are crossed. The change in sound speed from a 
value Clow to a value Chigh (and vice versa) is made on precisely the 
same time-scale as the change in viscosity parameter. Hence, no 
scaling is required for the amplitude of the luminosity variations 
driven by viscous dissipation in the disc. 



3.2 Simulations 

We perform simulations of an accretion disc in a binary stellar 
system with the parameters of the d warf n ova IP Pegasi taken 
from the catalogue of Ritter & Kol3 <1998h : Porb = 0.158 d. 
Ml = 1.15M0 and M2 = 0.67Mq. The accretion disc is built 
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Figure 1. V-band light curve for a simulation of the dwaii' nova IP Peg. 
Here, the inner part of the disc remains in the hot, high viscosity state 
throughout the quiescent intervals. The size of the hot region varies during 
quiescence, leading to the small variations in the light curve, but keeping the 
mean quiescent level flat. The four times t\ to are defined in the caption 
of Figure 2. 




Figure 2. Evolution of surface density in the accretion disc during quies- 
cence in the first simulation. The curves correspond to the following times 
in Figure [1] ti=26 (dots), i2=30 (dashes), (3=34 (long dashes) and t4=38 
(dot-dashes). The solid straight lines show the critical surface densities used 
in the simulation. 



up from scratch using a stream of particles injected from the in- 
ner Lagrangian point at a constant rate of 10^^" Moyr"^. We use 
Chot = 18.7 gcm~^, Ccoid = 7.02 gcm~^, Ttriggcr = 1750 s, 
viscosity parameters (^c = 1 and (^h = 10, corresponding to Oc = 
0.125 and Oh = 1.25, and sound speeds ciow = 0.05aQb and 
Chigti = 0.15anb. 

Figure Q shows the resulting V-band light curve. The light 
curve is calculated under the assumption that the disc is optically 
thick and each annulus of gas radiates as a black body. The Planck 
function is integrated over the interval 500 to 600 nm and summed 
over the annuli. In this simulation, many particles in the inner part 
of the disc remain constantly in the hot, high- viscosity state. In this 
region the critical surface densities given in the equations E| are 
both low and close together. The overall mean brightness is roughly 
constant during quiescence, with small variations occurring as the 
size of the hot region changes. It can be seen that this behaviour 
is now consistent with that observed in dwarf in quiescence. Our 
predicted variability in quiescence has too small an amplitude to 
be detected the visual light curves of amateur astronomers, but it is 
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Figure 3. Evolution of the contributions to the viscous energy dissipation rate made by the high and low viscosity particles in the disc for each of the times ti 
(top left) to it4 (lower right). The solid curve shows the contribution from low viscosity particles ^ 1.5), the dashed curve shows the contribution from high 
viscosity particles (f > 1.5), and the dotted curve the total viscous energy dissipation rate. The increasing low- viscosity contribution is offset by a decrease in 
the contribution from the high viscosity particles. 



Strikingly similar to that observed in the few detailed CCD observa- 
tions o f dwarf novae in quiescence (see, for example lWagner et alJ 

Figure |2| shows the variation in surface density during one of 
the quiescent intervals. It can be seen that the surface density in the 
inner disc never drops below Emin and hence remains in the hot 
state throughout quiescence. 

A hot inner disc seems to be a natural feature of these simula- 
tions, and in fact is a natural consequence of the linear triggers that 
get progressively closer together in the inner disc. With triggers of 
this form it is difficult to see how the inner regions of the disc can 
avoid entering the high state during quiescence. It is not immedi- 
ately obvious to us why this is not a feature seen in one-dimensional 
models. 

The response of the disc is shown in more detail in Figure|3] 
which shows the evolving contributions to the local viscous energy 
dissipation rate of the low and high viscosity particles at each ra- 
dius in the disc. For the purposes of this figure we define the high 
viscosity particles as those with C, > 1.5, to distinguish between 
those particles with a genuinely low viscosity, such as those arriv- 
ing from the stream which all have <^ = 1.0, and those which may 
be changing viscosity from a recent high state. At the beginning 
of quiescence, when the surface density and dissipation rate in the 
outer disc are low, the high viscosity particles dominate the inner 
disc. The hot particles at large radii [R > 0.2a) in panels 1 and 2 re- 
main from the previous outburst. By mid-quiescence (panel 3), they 



have cooled completely and the outer disc is dominated by the cold 
particles arriving from the mass stream. As the cold, outer regions 
begin to refill, the corresponding increase in dissipation is offset by 
a decrease in the size of the residual high viscosity inner region. An 
estimate for the transition radius, can be made from the inter- 
section of the solid and dashed curves in Figure|3| Although iih is 
varying both up and down during quiescence, the general trend is a 
decrease from i?h ~ 0.15a in panel 1, to i?h ~ 0.08a in panel 4. 

Throughout the quiescence, the outer radius shrinks slightly, 
in accordance with the expectations discussed in Section 2. We find 
that over the course of quiescence, although the mean luminosity 
remains constant, the total disc mass increases by 19 %, which is 
comparable to the mass accreted by the primary in each outburst. 

A second simulation of IP Peg was performed in order to 
demonstrate that it is the hot inner parts that keep the mean qui- 
escent brightness constant. In this case, we used critical surface 
densities that did not have a linear radial dependence, but were 
constant: Emax = Cup and Emin = Cdown, where we used 
Cup = 4.1 gcm^^ and Cdown = 1.6 gcm~^. In this case, it is 
more unlikely that the surface density of the disc in the inner parts 
can reach the critical level during quiescence, so the inner part of 
the disc remains in the cool, low-viscosity state. 

Figures |4] and |6| show the resulting V-band light curve and 
variations in surface density during the second simulation. The out- 
bursts are almost identical in profile, and the more familiar solution 
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Figure 4. In this simulation of outbursts in IP Peg, tlie suiface density trig- 
ger for the disc instability in the inner part of the disc is raised so that this 
region cannot reach the hot, high viscosity state during quiescence. Con- 
sequently, the mean brightness level rises steadily between outbursts. The 
four times ti to 44 are defined in the caption of Figure 6. 



is recovered, with a mean quiescent level that steadily increases 
between the outbursts as gas accumulates in the disc. 

The increase in V-magnitude during quiescence (~ 0.2 — 0.3 
mag.) is smaller than the increase seen in 1-d models 1 mag.). 
This can be understood in terms of the sound speed, which only 
changes (that is, the disc only heats up) when the viscosity changes. 
In the 1-d models, the temperature and sound speed of the disc 
increase during quiescence as the surface density increases, even 
when a = ftcoid everywhere. It is trivial to estimate the effect that 
this change in sound-speed would have on our simulation. The dis- 
sipation rate, D{r) ~ !/E ~ TcE, where Tc is the mid-plane tem- 
perature. In our simulation, the dissipation rate at the end of quies- 
cence is about 1.25 times (0.2 magnitudes larger than) the rate at 
the beginning of quiescence. This arises from the increase in sur- 
face density alone. Given a reasonable increase of a factor of two 
in Tc over quiescence, the increase in dissipation rate would be a 
factor of two higher, that is a difference of 2.5 times (or 1 magni- 
tude). Indeed, we can go one step further and incorporate this into 
the simulations. 

In an a-disc, surface density scales with sound speed as E ~ 



3/14 



For all particles with a = Qc, we let 



Cs — Clo 



(16) 



The results can be seen in Figure|5| For the flat-trigger case (sim- 
ulation 2), where most of the particles during quiescence have 
Q = Qcoid, the rise in the V-band light curve is now a full mag- 
nitude, in line with our expectations. The difference in simulation 1 
is much less marked, because there the viscosity is not Ocoid every- 
where during quiescence. The amplitude of the variations in quies- 
cence in the first simulation remain much smaller than the increase 
in V-band magnitude in the second simulation. 

Figure Q shows the evolution of total viscous energy genera- 
tion over the disc in the second simulation. Again, the outer radius 
moves inward slightly, but there is a general increase in dissipation 
rate at all radii. The increase in disc mass is virtually identical to 
that found in simulation 1, around 19 %. 
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Figure 5. Evolution of the V-band light curve when the sound speed is var- 
ied with surface density for the cold particles. Top: the light curve of the 
second (flat-trigger) simulation now shows a brightening by a full magni- 
tude in quiescence. Bottom: the difference in the first simulation is much 
less marked. 
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Figure 6. Evolution of surface density in the accretion disc during quies- 
cence in the second simulation. The curves correspond to the following 
times in Figure|l] ti=24 (dots), t2=29 (dashes), 1^=34- (long dashes) and 
t4=39 (dot-dashes). The solid straight lines show the critical surface densi- 
ties used in the simulation. 



4 DISCUSSION 

We have shown that the observed constant mean brightness of 
dwarf novae during quiescence can be regulated by the response of 
a small inner portion of the disc that remains in a high-viscosity 
state close to the critical surface density limits for the thermal- 
viscous disc instability. Over the course of quiescence, the increase 
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Figure 7. Evolution of total viscous energy dissipation rate in the second 
simulation. The lines correspond to the same times as those in Figure 6. 



dwarf novae, there is observational evidence to support the pres- 
ence of such a hot, optically thin inner region (see , for example, 
llaDous et al]^1997^ : lDone C. & Osborne J.Pj<1997h l Such a flow 
would almost certainly be highly turbulent, and it is exactly these 
conditions under which the MRl is expected to be most efficient 
and the viscous stress high. 

We conclude that the constant brightness level that is observed 
between the outbursts of dwarf novae is well explained by the grad- 
ual cooling of a small, critically-stable, hot inner region. The di- 
minishing contribution to the viscous energy dissipation rate from 
such a region is sufficient to mask the increasing energy dissipation 
due to increases in mid-plane temperature and surface density in 
the rest of the cool disc. 
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in dissipation in the outer disc brought about by the increasing sur- 
face density and central temperature is offset by the gradual cooling 
and shrinking of the high-viscosity inner region. 

A steady increase in the mean brightness level between 
outbursts is unavoidable if the entire disc is in a single vis- 
cosity state, because both the surface density and temperature 
must be increasing. This has long been known to be the case 
from thermodynamically-detailed one-dimensional calculations in 
which the quiescent disc is entirely in the cool, low viscosity state. 
The effect of increasing E and T dominates the effect of the 
disc shrinking in quiescence (equation^, and the brightness rises 
steadily. 

However, we can see from eg nation fTol that this problem can 
be resolved if a small region remains in the high-viscosity state 
while the surface density and temperature increase in the rest of 
the disc. In the SPH calculations, this is achieved by the surface 
density hovering near the critical stability limit at all times. The 
high-viscosity, hot inner disc is the main difference between our 
simulations and t hose within th e framework of the standard disc 
instability model Casot^'^OoT), where the quiescent disc is cold 
everywhere and the surface density is left well below Emin after 
each outburst. 

The key question that will remain is how can this be achieved 
in a real accretion disc? It has been suggested that the inner parts 
of the ac cretion disc can be heated and evaporated near the white 
dwarf tMever & Mever-Hofmeister 1994,) . Indeed, a low density 
region in the inner disc can explain the UV delay observed in dwarf 
novae, because at the onset of an outburst it takes a few hours for the 
low density region to refill. Of course, the UV delay will be present 
however the low de nsity region arises: in addition to the evapora- 
tion m odel of . Meyer & Mever-Hofmeisteilh994) . lLivio & Pringld 
il992h have suggested that the inner part of a quiescent disc can 
be evacuated by the weak magnetic field of the white dwarf pri- 
mary. An attractive feature of the models that we have presented 
in this paper is that a similar low-density inner disc is predicted 
during quiescence. Here, however, the small low density region is 
a consequence of a local viscosity higher than that which is usu- 
ally assumed for a quiescent disc. This is exactly analogous to the 
situation found by Stehle & King ( 1998), who showed that the UV 
delay can be reproduced if irradiation by the hot white dwarf keeps 
the disc hot and fully ionised out to several stellar radii. 

Although there are difficulties separating spectral components 
due to the accretion flow itself and the white dwarf in quiescent 
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